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1 Introduction 



As is well known, all structure functions of the lowest twist-two, F\(x) , F 2 {x) , g\(x) which are 
measured in the deep-inelastic lepton-nucleon scattering, conserve chirality. Ralston and Soper 
[|IJ first demonstrated that besides these structure functions, there exists the twist-two chirality 
violating nucleon structure function h\{x). This structure function does not manifest itself in 
the deep inelastic lepton-hadron scattering, but can be measured in the Drell-Yan process with 
both beam and target transversely polarized. The reason of this circumstance is the following. 
The cross section of the deep inelastic electron(muon)-hadron scattering is proportional to the 
imaginary part of the forward virtual photon-hadron scattering amplitude. At high photon 
virtuality the quark Compton amplitude dominates, where the photon is absorbed and emitted 
by the same quark (Fig. la) and the conservation of chirality is evident. The cross section of 
the Drell-Yan process can be represented as an imaginary part of the Fig. lb diagram. Here 
virtual photons interact with different quarks and it is possible, as is shown in Fig. lb, that 
their chiralities are also different. It is clear from Fig. lb, that chirality violating amplitude 
in Drell-Yan processes is not suppressed at high Q 2 in comparison with chirality conserving 
ones and, consequently, corresponds to twist two. However, this amplitude corresponding to 
target spin flip, has no parton interpretation in the standard helicity basis. Nevertheless, h\{x) 
can be interpreted as a difference of quark densities with eigenvalues +1/2 and —1/2 of the 
transverse Pauli-Lubanski operator Sj_75 in the transversely polarized proton as it is explained 
in details by Jaffe and Ji in ||. 

Until now, there is no experimental data on the chirality violating nucleon structure func- 
tion h\(x). Besides |J1 0] the theoretical study of this structure function ( under the name 
Aiq(x) ) have been performed also by Artru and Mekhfi ||. The role of h\(x) in factorization 
of a general hard process with polarized particles was investigated by Collins M. The first 
attempt to calculate h±(x) was carried out by Jaffe and Ji by means of the bag model. 

In this paper we calculate the proton structure function h\{x) by means of the QCD sum 
rule approach. The idea of the method is the following [||, |6j. Consider the four-point vacuum 
correlator corresponding to the amplitude of forward scattering of the current r)(x) with proton 
quantum numbers on external currents j\ (x) , j 2 (x) . (In the case of the structure function 
F 2 (x) considered earlier |^|, || ji(x), j 2 (x) were electromagnetic or weak currents). Let the 
momenta, corresponding to the currents f]{x) and j(x) be p and q. As was shown in [^J, |6|, 
the imaginary part of the forward scattering amplitude is determined by small distances in 
t-channel if p 2 and q 2 are negative and large enough, | p 2 |, | q 2 |^> R~ 2 , (R is the confinement 
radius) and the Bjorken scaling variable x = Q 2 /2z/, Q 2 = — q 2 , v = p ■ q is not close to 
the boundary values x = and x — 1. Therefore, to calculate the amplitude of interest we 
may use the operator product expansion (OPE) method accounting the vacuum expectation 
values (v.e.v.) of various operators. If, in addition, we suppose | p 2 |<C| q 2 | and restrict 
ourselves to the first term in expansion over p 2 /q 2 , only twist-two structures will be retained. 
In order to obtain the desired proton structure function we use the standard procedure of the 
QCD sum rule approach. The amplitude calculated in QCD is equated to the contribution of 
physical states, the proton and the excited states with proton quantum numbers, by means of 
the dispersion relation in p 2 . The proton contribution we are interested in, is then separated 
by applying the Borel transformation in p 2 which suppresses the contribution of excited states. 

It is important to check that the results satisfy necessary conditions of QCD sum rule 
calculations: 1) the OPE is converging i.e., the the highest accounted term in OPE is well 
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less than the first term(s) of OPE and 2) in the phenomenological part of the sum rule the 
contribution of excited states does not surpass the one of the proton. These requirements 
determine the domain of x, where the results are valid. The value of Q 2 must be much larger 
than the proton mass m squared, Q 2 3> m 2 , but cannot be very large, since the leading 
logarithmic corrections in QCD are not accounted in our approach. So, we expect that our 
results are valid at Q 2 ~ 5 — lOGeV 2 . 

It should be emphasized that calculation of various proton structure functions by means 
of QCD sum rules contains no experimental inputs. Moreover, there are no new parameters 
at all. The values of vacuum condensates determining the OPE are fixed from QCD sum rules 
for various two-point correlators, whereas the parameters of the proton channel, namely the 
transition constant of proton into the quark current and the value of the continuum threshold 
are taken from the analysis of the sum rules for the proton mass 0. The price for model 
independence is a limited domain of x where the calculated structure function is reliable. 
Thus, for the valence u(x) distribution calculated in || the region of validity 0.2 < x < 0.7 
was established. Calculated by the same method, the u quark contribution to the proton spin 
structure function g\{x) || has an even narrower interval of validity 0.5 < x < 0.7. The 
criteria determining these intervals as well as the accuracy of the calculation are well defined 
from comparing nonleading and leading terms of OPE, as explained above. A comparison 
with experiment shows that one can rely on this method of calculation within the domain of 
validity. A very nontrivial evidence of that is the accordance of the calculated value of g\{x) 
with recent experimental data ( see detailed discussion below). 

In carrying out this program for calculation of h\(x) some specifics arise. In refs.JT], [| 
the structure function h\(x) was represented in terms of the matrix elements of the operators 
on the light-cone. In our calculation we start from the forward scattering amplitude on the 
proton. In Sec. 2 the forward scattering amplitude corresponding to the light cone definition 
of hi(x) is constructed. The four-point current correlator calculated in QCD contains many 
spin-tensor structures and it is necessary to find such structure, which, if saturated by the 
proton state, gives the desired function hi(x). This is done in Sec. 3. 

In sec. 4 we calculate the OPE of the four-point current correlator. Here the most important 
difference between the function h\(x) and other twist-two structure functions appears: since 
h\ (x) violates chirality, the OPE in the chosen spin-tensor structure starts from the operator 
of dimension 3 - the quark condensate, unlike the case of other twist-two structure functions, 
where it starts from the zero dimension unit operator. In addition to the main term of OPE 
we calculate the next term proportional to the mixed quark-gluon condensate of dimension 5. 

In Sec. 5 we obtain the sum rule, perform numerical analysis and give the results of the 
hi(x) determination. Our conclusions are presented in Sec. 6. 

2 The connection of light-cone and forward scattering 
amplitude definitions of the structure function h\(x) 

The structure function h\ (x) was defined in ref. |2[] by Fourier transformation of the operator 
product on the light cone between proton states 
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^e iXx <p,s\ ^(0)<v7 5 V>(An) \p,s>= 

2[hi{x, q 2 ){s^p u - sj_ u p^) + h L (x,q 2 )m 2 (p^n u - 

p v n^)(s ■ n) + h 3 (x, q 2 )m 2 (s±^n u - s± u n fl )] (1) 

Here n is a light cone vector of dimension of (mass) -1 , n 2 = 0, n + = 0, p-n = 1, p and s are the 
proton momentum and spin vectors: p 2 = m 2 , s 2 = —1, p-s = 0, and s = (s-n)p+(s-p)n+s±. 
We choose in this section the proton and virtual photon c.m.system and put q = {qo, 0, 0, E}, 
p = {E, 0, 0, —E}, n = {1/E, 0, 0, l/E}/2, where E is the virtual photon momentum in this 
system and the terms of order m 2 /E 2 are neglected. 

Let us demonstrate now that h\(x) can also be defined as a matrix element between proton 
states for forward scattering of the axial current with its transition into scalar current pA —>■ pS 
plus crossing matrix element 

T^(p,q,8, 8 ) = % J d 4 xe iqx <p,s\ (1/2)T{j m5 (z),.7(0) +j(x)J^(0)} \p,s> (2) 

where j^(x) and j(x) are axial and scalar currents. Supposing approximate conservation of 
axial current, T^(p, q, s, s) can be represented as 

q ' s ~ 

r M (p, q, s, s) = (s M - —q^hix, q 2 )+ (3) 

+ - —jrKQ ' s )h(x, q 2 ) + e^xaPuqxSaiq ■ s)l 2 (x, q 2 ) 
(only spin-dependent terms are retained). We will prove that 

h\(x, q 2 ) = Imhi(x,q 2 ) (4) 

7T 

up to higher twist terms. 

Until now we did not specify the flavors of axial and scalar currents. Strictly speaking, 
due to the anomaly, the representation (3) is correct only for flavor octet axial current. For 
the singlet axial current additional terms appear in (3). But, as we shall see, in proton the 
contribution of w-quarks to hi(x) is much larger than of ci-quarks. Therefore, there is no 
matter whether we consider uu current, as we shall do, or uu — dd current, where anomaly 
does not contribute and (3) is correct up to the terms of order m 2 /Q 2 . In any case, the 
anomaly contribution is suppressed by a s , what is beyond the accuracy of our results. For 
these reasons, we use eq.(l) for the case of one flavor. 

As is well known, at high Q 2 the imaginary part of the forward scattering amplitude like (2) 
is dominated by the light cone region in space-time. Using the expression for quark propagator 
near light cone 

< | T{^(x),^(0)} | >= --^( 7 ^A_L_ ( 5 ) 



and substituting it into (2) we obtain 



ImTf, = -— I d 4 xe iqx <p,s\ V>0)(z7m - 7^)75^(0) \p,s> 5'{x 2 ) ( 6 ) 



1 We are indebted to R.Jaffe for the idea of this definition of hi(x). 
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Using the light cone variables 

t-z = r (l/2)(t + z) = X (7) 

we can rewrite (6) as 

ImT^ = - l -Jd\JdrJ dx 2 e lEiT - 2 ^6(2T\ - x 2 )6'(x 2 )x (8) 

x < p, s | i)(x)a tlv x v j 5 '4)(0) \p,s> 

where the notation £ was temporarily accepted for the Bjorken scaling variable and we used 
the relation 

Q 2 

qx » E(t -z) + £=t « Er- 2£EX (9) 
2E 

Integrating over x 2 and r in (8) and putting on the light cone 

x u = 2E\n v (10) 

we get 

oo 

ImT^ = ~~J d\'e iX 'Z <p,s\ ^)a^n v ^{nX') \p,s>, (11) 

— oo 

where A' = 2EX. Equation (11) can be compared with (1), multiplied by n u . The comparison 
of terms, proportional to s M , gives eq.(4). It should be mentioned that instead of axial and 
scalar currents one can also use a combination of vector and pseudoscalar currents j At (x)j5(0) — 



3 Selection of spin-tensor structures 

Consider the four-point correlator 

n„(p,g) = -ij d 4 xd 4 yd 4 ze iqx+ip{y ~ z) x (12) 

x < | T{ V (y),l/2(j^(x),j(0) +j(x),j tl5 (0)),fj(z)} | > 
where 77 is the three-quark current with the proton quantum numbers M 

V = e abc (u a C lx u b ) l57x d c (13) 

a,b,c = 1,2,3 are colour indices. We are interested in the proton contribution to U^(p,q), 
which is given by 

UP(p,q) = X% 1 Yl v r (p)TW(p, q, r, r')v r ' (p) (14) 

\P m ) r y 

where v r (p) is the proton spinor with momentum and polarization r, A ^ is the transition 
constant of proton into quark current 
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< I i] I p, r >= \nv v (p) 

and Tjp(p, q,r,r') is nondiagonal in the proton spin matrix element (2). It is convenient to 
represent Tjf\p, q,r,r') in the form 

Tjf\p,q,r,r') = v r (p)f^(p,q)v r '(p), (15) 

Here Tjf ^ (p, q) is the matrix element between proton states before multiplication on proton 
spinors, i.e. where p 2 = m 2 , but p ^ m. 

Since the axial and scalar currents enter symmetrically (2) and (12), Il M (j9, q) as well as 
Tjf^ (p, q) satisfy the condition 

F(q,p)=C' l F T (-q,-p)C, F = U„ (16) 

where C is the charge conjugation matrix. Accounting for (16) the general structure of 
U'Hp.q) is 

f ^(p, q) = [K{q ltl - 1 ,q) + E(pq - qp) Ptl + Cp,q + C'q.p (17) 

+D(pq - qp)q^ + L(p7 M - 7^) + M^^ + Rqq^ + Spp^ 

where K, E, C, C, D, L, M, R, S are functions of v = p ■ q, q 2 and p 2 = m 2 . n M (p, q) has the 
same structure, in this case p 2 being not equal to m 2 . 

The matrix element of Tjf^ (p, q) over the polarized proton state with spin s M is given by 

(l/2)Tr{f ^(p, q){\ - 7s s)(p + m)} = -As^uK + Lm 2 + (1/2) Mm] + ... (18) 

where only the terms proportional to s M are written explicitly. In accord with (3) these terms 
are proportional to hi(x). Substitute (15) into (14) and use the relation 



We have 



r 

n? = { /} m ,/ P + ™)TP(p, q)(p + m)= ^ ^ 2)2 x 
A/ff) 

{ — —[{yK + Lm 2 + (l/2)Mm) + m 2 q 2 (D + (l/2)R)\p^q 

+ 4m 3 (D + (1/2)^)^)75 + ••• (19) 

where the terms at the structures p^q^ and q^qj^ in the decomposition (17) are separated. 
When deriving (19) we used the relation IT M g M = following from conservation of axial current. 
As is seen from (18), (19) the combination of the structure functions vK + Lm 2 + Mm/2, 
which is proportional to hi (x) can be found by studying U.^ (p, q) as a coefficient function 
at the structure p^q^s added by the coefficient function at the structure q^q^, multiplied by 
q 2 jv = —2x. In fact, the latter, as can be easily seen, and will be proved by direct calculation 
in Sec. 4 corresponds to twist 4 and can be neglected. 

Therefore, the recipe of the calculation of h\(x) is the following. Calculate 
JmlLjQo, q) defined by eq.(12) in QCD at p 2 < using OPE in 1/p 2 and retaining the leading 
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terms in 1/q 2 , corresponding to twist 2 amplitude. Separate the coefficient function at the 
structure p^p'js- On the other side, using the dispersion relation in p 2 represent the same func- 
tion through the contribution of the physical states - the proton and excited states, the latter 
approximated by continuum. Equate two representations and apply the Borel transformation 
to both sides of the sum rule in order to enhance the proton contribution and improve the 
convergence of OPE. This (up to some details) gives the desired sum rule for determination 
of the proton structure function hi{x). 



4 QCD calculation of four-point correlator 



We calculate the imaginary part in s-channel of the four-point correlator n M (p, q) (12). Since 
H^(p,q) violates chirality, the OPE starts from the term proportional to quark condensate 
< | qq | >. The corresponding diagrams are given in Fig. 2. 

Consider first the case, when the axial and scalar currents are w-quark currents (Fig. 2a). 
The calculation of the diagram Fig. 2a gives 

ImUVp, q) = --— < I mm I > f(-4 + 9x+ 

2 

— p X 

+ 2xln^—)p fI q + g M g]7 5 (20) 

(only the structures p^q^s and q^q are retained). After the Borel transformation in p 2 the 
contribution of the structure q^q^ vanishes and we get 

B M 2lmU u Jp, q) = < | uu | > M 2 xp^q l5 (21) 
Lixv 

where M 2 is the Borel parameter. 

In the case, when the scattering proceeds on (i-quark, it is clear that the corresponding 
amplitude (Fig.2b) in the lowest twist depends only on v and q 2 and is p 2 independent. (Direct 
calculation gives the value — (l/127r) < | dd | > jvx for the term proportional to p^cfi^)- 
Thus, it vanishes after Borel transformation. The nonvanishing contribution in the case of 
scattering comes from diagrams with additional gluon exchange in the diagram of Fig. lb. 
They are of the order (a s /ir) times the r.h.s. of eq.(21) and we will disregard them, since they 
are of the same order as perturbative corrections to eq.(21), which are neglected. 

The next term in the OPE for H^(p, q) is proportional to the mixed quark-gluon condensate 
of dimension 5: 

- g < | qa^G^q | >= m 2 < | qq | >, (22) 

Consider again first the case, when the scattering proceeds on w-quarks. The calculations 
will be performed in the fixed point gauge 

x,A;(x) = (23) 

and we choose as a fixed point the left-side lower vertex point in Figs. 1,2. The first source of 
appearance of quark-gluon condensate is the expansion of quark condensate in the diagram of 
Fig. 2a in powers of x. We write now 
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< I T{<(x), 4(0)} | >= -(l/12)5 ab 8 af3 < | u(0)u(0) | > + 

+ Y^ 5ab5 <xpg < I u(0)°^\ n u(0) I > x 2 (24) 

The first term in the r.h.s. of eq.(24) gives the already accounted quark condensate con- 
tribution. After substituting the second term in (24) into the diagram of Fig. 2a we obtain for 
the coefficient function at the structure p^q^s 

ImU^ = §■ < | | > -2V375 (25) 

VKV p z 2 

The structure g^^Ts is absent in twist 2 terms. In other way the quark-gluon condensate 
emerges from the diagrams of Fig. 3. To calculate these diagrams we used the quark propagator 
in the constant gluonic field ||: 
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+ -gX-G n a ^( la k 2 -2k a k)-} (26) 

The results of the calculations of the Fig. 3 diagrams are the following. The diagrams with gluon 
emission from both horizontal lines of u and <i-quarks (Figs.3a,b) give a vanishing contribution 
to the coefficients at the structures we are interested in. The sum of the diagrams Figs. 3c, d 
with gluon emission from vertical w-quark lines, when the second term in the r.h.s. of (26) is 
accounted, gives 

I TTV^ 1 1 — CC 

ImU^ = — £ < | uu | > ~( )^g 75 (27) 

Attu p z 6 x 

The most complicated is the calculation of the third term in the r.h.s. of the (26) contribution. 
In the chosen fixed point this term contributes only to the Fig. 3d diagram and the calculation 
gives: 

i tti 2 i 

= -- £ < | im | > — Pfl q"f 5 (28) 

Attu p z 12 

For all Fig. 3 diagrams the coefficient function at the structure (^75 vanishes in twist 2 
terms. The total contribution of dimension 5 mixed quark-gluon condensate to the four-point 
correlator, when the scattering proceeds on w-quarks, is obtained by summing (25), (27), (28) 

1 m 2 11 1 

Iml% = §■ < I uu I > -(- - - - 3x)p^q l5 (29) 

p Airiy p A bx2 

After borelization we have 

B M 2lmU u (p, q) = -^—ml <Q \ uu\Q >\(- - \ - 3x) Pfl q-f 5 (30) 
p Airu 6 x 2 

For the case of scattering on <i-quark the contribution of quark-gluon condensate vanishes for 

the same reason as the quark condensate contribution. Since in the scattering on <i-quark 

the first two terms in OPE vanish and the charge square of <i-quark is 4 times smaller than 
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w-quark, we can safely disregard the <i-quark contribution to the proton structure function 
ht(x). 



5 The sum rule. Determination of h\(x). 

In accordance with the results of Sec. 3 consider the coefficient function H(p 2 ,x) at the struc- 
ture p M g75 in the four-point correlator (12) as a function of p 2 . Using the double dispersion 
relation represent H(p 2 ,x) in terms of the contribution of physical states - the proton and 
excited states 

tt( 2 \ \ 2 N inn A(x) °f p(pl,x) 2 

H(p,x) = — —h 1 (x) + - + / — —d Pl 31 

v (p 2 — m 2 ) 2 v p 2 — m 2 J 2 {pi — p 2 ) 2 

The factor — it in the first term in (31) comes from (4) and the factor mjv emerges from 
comparison of (18) and (19). The last term in (31) corresponds to the contribution of excited 
states approximated by continuum. In our case this contribution is determined by the diagram 
of Fig. la and is given by eq.(20), W 2 is the continuum threshold. The second term in (31) 
is a nondiagonal term where the current f] produces a proton from the vacuum, which, after 
interaction with and j goes into some excited state N* absorbed finally by the current 
n (for details see ||). Apply the Borel transformation to (31) and equate it to the result of 
QCD calculations (21), (30). We obtain the sum rule for w-quark contribution to the proton 
structure function hi(x) 

Uu( \ m 2 /M 2 ti/2 

W +B <^ 2 %^ 2M2 ^<^ 4/9 < 32 > 



where B = A/X 



2 

N ■ 



a = -(2tt) 2 < | qq | >= 0.55GeV 2 
X 2 N = 32tt 4 X 2 n = 2.1GeV G , W 2 = 2.3GeV 2 (33) 
E x {z) = l-e- z (l + z) 

and powers the of factor L = ln(M/A)/ln(fi/A) take into account anomalous dimensions 
of currents in the correlator and operators in the OPE. ( see ]F], for details). Note that 
the anomalous dimension of the scalar current cancels with the corresponding anomalous 
dimension of the quark condensate operator in the first term of (33), whereas the anomalous 
dimension of the quark- gluon condensate operator is small and we neglect it. We assume A = 
150MeV and fi = 0.5GeV. The numerical values in (33) correspond to the best description of 
nucleon mass and magnetic moments in QCD sum rule approach ||, [TDf] . In order to get rid 
of unknown contribution of B(x) let us differentiate (32) over 1/M 2 . We get 

K(x) = 2^e m2 / M2 {2M 2 x[(m 2 - M 2 ) El (^) + ^ e -^/M 2 ]L -4/ 9+ 
mX 2 N M 2 M 2 
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+ I m 2 m 2 ( I_ 1 _ 3x ) L -8/9 } (34) 

As is seen from (32), (34), these expressions are not correct at x — > 0, the OPE breaks down 
at small x, as is expected from general grounds f6j. In order to find the values of x, where we 
can believe our results, compare the magnitude of the second term of OPE with the first one. 
We find that at x > 0.4 the second term of OPE comprises less than 20-25% of the total in 
eq.(34) and less than 15-20% in eq.(32). Therefore we may have a confidence in our results at 
x > 0.3. At x close to 1 our results are not reliable because it is a resonance region (see ref. 
[BJ for more detailed explanation). In Fig.4 we plot the Borel parameter M 2 dependence of 
/i"(x) determined according to eq.(34). As is seen from the Figure, the dependence is rather 
weak and this necessary condition for the QCD sum rule validity is fulfilled. Our basic result 
- the M-quark contribution to proton structure function hi(x) determined according to (34), 
is plotted solid curve in Fig. 5 (at M 2 = lGeV 2 ). 

It must be mentioned that differentiation over 1/M 2 , which was used in order to get (34) 
from (32), spoils the accuracy of the sum rule, since the role of continuum and of unaccounted 
higher order terms of OPE increases. Particularly, the continuum contribution, which com- 
prises 30% in eq.32 is about 70% or even more, especially at large x in eq.(34). (these values 
refer to M 2 = lGeV 2 ). For this reason and also because we accounted only two terms of OPE, 
the accuracy of our results is not high. We estimate it by ~ 30% at x ~ 0.5 and by 50% at the 
border of interval where our approach is working, x = 0.3 and x = 0.7. At a; > 0.7 the curve 
of h^(x) given in Fig.4 is not correct, since it violates the inequality | h\(x) \< u(x) proved in 
ref. |. 

As was explained above, the contribution of <i-quarks is small and the proton structure 
function 

h\{x) « (4/9K(x). (35) 

In Fig. 5 we also plotted h±(x) found by Jaffe and Ji in the bag model [Q] (dashed curve). 
The bag model curve is 2-3 times higher compared with the results of our calculations and 
strongly depends on x at intermediate x, while our hi(x) is flat (due to errors we cannot insist 
on negative curvature of h^(x), probably it is an artifact of our approximations). It should be 
mentioned that the bag model curve strongly violates the inequality hi(x) < u(x), if one uses 
the experimental values for u{x). 

In order to establish an input for comparison with data of future experimental studies as 
well as for possible renormalization group analysis we suggest the following ansatz for h\{x) 
in the whole region of Bjorken variable, based on our calculation. We recall first that the 
behaviour of structure functions at x — > is governed by the Regge trajectory in the t-channel 



of the corresponding forward scattering amplitude ( see for details, JTT]]) so that, in particular, 

h{(x) ~x~ a ^ 0) (36) 

where a ai (0) is the intercept of the axial ai-meson trajectory. Assuming that this trajectory 
is linear and has the same slope as p, a2 trajectories, namely a' ~ !GeV~ 2 , and using the 
mass of the ai meson we easily get a ai (0) ~ —0.3. We then extrapolate h\(x) from the point 
x = 0.3, the minimal value at which our calculation is reliable, to the region < x < 0.3. 
using our result for h±(x = 0.3) and (35), (36). Note that Regge behaviour was also invoked 
in for analysis of the moments of h\{x). 
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In the large x region where the structure function should decrease according to the general 
quark counting rules, we simply use the inequality u(x) > h^(x) established in and replace 
h x (x) by the corresponding values of the u(x) distribution ( taking it, e.g. from |12| )beginning 
from the point x = 0.55 where h x (x) starts to violate this inequality. ( Within the uncertainties 
of h x (x) and u(x) this point may be as large as x ~ 0.7, the upper limit of the interval of 
reliability of our calculations). The resulting curve is shown in Fig. 6. 

Finally, let us compare our predictions for hi(x) with the nucleon spin structure function 
gi(x). Using recent data of SLAC E142 experiment [13| on the neutron spin structure function 



g™(x) measured at Q 2 > lGeV 2 we conclude that at x > 0.3 this structure function values are 
consistent with zero. In terms of separate u and d quark spin distributions in the proton, the 
following relation should be valid starting from x = 0.3 : g x (x) — —Agf(x) clearly indicating 
that the the proton spin distribution gi(x), as well as h\(x), is dominated at x > 0.3 by the 
u-quark contribution: 

g\{x) ~ (4/9)tf (x) . (37) 



Taking experimental results on gi(x) from EMC measurement [14] we get, e.g. g±(x 



0.5 (0.2) at x = 0.3 (x = 0.5). The recent SMC measurement QT5j gives the values of g x (x 



at x = 0.3 and 0.5 about 20% lower, but agreeing with EMC within errors. Comparison with 
our calculation clearly indicates the validity of the inequality for u quark distributions in the 
proton: 

h u x {x) > g u x {x) (38) 

at x > 0.3. The same conclusion can be obtained, if we compare h\(x) with the structure 
function g\{x) calculated in [Q in the framework of the method, similar to used here. Unfor- 
tunately, the gi(x) calculation is legitimate only in the narrow domain of x, 0.5 < x < 0.7. 
In this domain the results of || are in good agreement with recent data of SMC . In the 



interval 0.4 < x < 0.7 the result of this experiment for the mean value of the proton spin 
structure function is gf(x) = 0.08 ± 0.02 ± 0.01 (see Table 1 in |15|). The prediction of ref. || 
for the same interval of x was g\{x) = 0.06 with an estimated error about 30% (some extrap- 
olation of the g\{x) curve from x = 0.5 to x = 0.4 is done). As was shown in [§], the u-quark 
contribution to the spin structure function dominates in this region of x, in accordance with 
d37|) . Using the results of || we can compare theoretical predictions of h x {x) and g x (x) at 
the point x = 0.5 (but not at x = 0.3). We find from || g x (0.5) ~ 0.14 which agrees within 
theoretical and experimental errors with the value obtained above from EMC experiment and 
with inequality (|3"8"D. The results of bag model calculations M also obey 



6 Conclusion 



In this paper a new method of calculation of the proton chirality violating structure function 
hi(x) at intermediate x has been developed based on the QCD sum rule approach. The 
obtained estimates for this structure function are valid at 0.3 < x < 0.7. Since we disregard 
the perturbative QCD correction, our results are valid at intermediate Q 2 ~ 5 — lOGeV 2 , 
where deviations from scaling are inessential. In our approach h\(x) is proportional to the 
quark condensate. We have calculated only two terms in the OPE. Nevertheless, once the 
method is established, the accuracy may be improved in future taking into account higher 
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order terms of OPE. We estimate the current uncertainty of our calculation at the level of 
30-50%. Therefore, having in mind that the similar calculation of the proton spin structure 
function g%(x) has reasonable agreement with experiment we consider our estimate of hF^x) 
as a reliable and useful guideline for experiments aimed at measuring this structure function. 

We found that h\{x) is rather flat in the interval 0.3 < x < 0.7 and not too much smaller 
than an upper limit | K{{x) \< u(x) derived in f2|. This circumstance gives a good chance for its 
experimental study. Another interesting feature of h\{x) is that the main contribution comes 
from w-quarks, the d-quark contribution being suppressed, besides of the d-quark charge, by 
a factor of order a s /7r ~ 0.1 or is given by the contribution of higher dimension operators 
in OPE, which, probably, are also small. For this reason we expect that for the neutron at 
intermediate x 

til(x) « (l/4)/i?(x) 

The comparison of our h\{x) with the same function calculated by Jaffe and Ji in the bag 
model [T2J shows that our hf(x) is by a factor of 2-3 smaller in the interval 0.3 < x < 0.5 and 
matches with h\{x) at x — 0.6. The condition | hx(x) \>\ gi(x) | found in the bag model 
is also fulfilled in our calculation. 
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THE FIGURES 
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Figure 1: (a) Deep inelastic lepton-hadron scattering, the chirality of quarks is conserved. 
Solid lines are quarks, wavy lines are virtual photons, R(L) denote right (left) chirality of 
quarks; (b) Drell-Yan process with chirality of quarks flipped. 
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Figure 2: Diagrams corresponding to quark condensate contribution to the four-point cor- 
relator ri M defined in (12): a) leading quark condensate contribution from the scattering on 
u-quark ;b) nonleading contribution from the scattering on J-quark. Wavy lines denote axial 
or scalar currents, external momenta are shown by arrows, lines with a cross denote the quark 
vacuum fields. 
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Figure 3: Diagrams corresponding to quark-gluon condensate contribution to the correlator 
(12). Dashed lines denote the vacuum gluon held. Black squares denote the hxed point of the 
gauge. 
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Figure 4: Dependence of u-quark contribution to hi(x) on Borel parameter M 2 at various 
values of x 
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Figure 5: The u-quark contribution to the proton structure function hi(x) determined from 
QCD sum rule (34) at M 2 = lGeV 2 ( solid curve). The bag model prediction obtained by 
Jaffe and Ji [2] for the same structure function is shown by dashed line. 
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Figure 6: The u-quark contribution to the proton structure function hi(x) based on QCD 
calculation at intermediate values of x. In the region of low x an extrapolation from the point 
x=0.3 is done according to the expected Regge behaviour of this structure function. The 
dashed line shows the valence u-quark distribution taken from ref. [12] at Q 2 = 5GeV 2 which 
serves as an upper limit for h™(x) at large x. 
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